lAoB Plants 



The open-access 
journal for plant sciences 



Invited Mini-Review 

SPECIAL ISSUE: Scaling Effects Regulating Plant Response 
to Global Change 



Plant community feedbacks and long-term ecosystem 
responses to multi-factored global change 

J. Adam Langley 1 * and Bruce A. Hungate 2 

1 Department of Biology, Villanova University, Villanova, PA 19085, USA 

2 Center for Ecosystem Science and Society, Deparment of Biological Sciences, Northern Arizona University, Flagstaff, AZ 86011, USA 
Received: 1 April 2014; Accepted: 18 June 2014; Published: 14 July 2014 

Guest Editor: Elise S. Gornish 

Citation: Langley JA, Hungate BA. 2014. Plant community feedbacks and long-term ecosystem responses to multi-factored global 
change. AoB PLANTS 6: plu035; doi:10.1093/aobpla/plu035 



Abstract. While short-term plant responses to global change are driven by physiological mechanisms, which are 
represented relatively well by models, long-term ecosystem responses to global change may be determined by shifts 
in plant community structure resulting from other ecological phenomena such as interspecific interactions, which are 
represented poorly by models. In single-factor scenarios, plant communities often adjust to increase ecosystem 
response to that factor. For instance, some early global change experiments showed that elevated C0 2 favours plants 
that respond strongly to elevated C0 2 , generally amplifying the response of ecosystem productivity to elevated C0 2 , 
a positive community feedback. However, most ecosystems are subject to multiple drivers of change, which can 
complicate the community feedback effect in ways that are more difficult to generalize. Recent studies have shown 
that (i) shifts in plant community structure cannot be reliably predicted from short-term plant physiological response 
to global change and (ii) that the ecosystem response to multi-factored change is commonly less than the sum of its 
parts. Here, we survey results from long-term field manipulations to examine the role community shifts may play in 
explaining these common findings. We use a simple model to examine the potential importance of community shifts 
in governing ecosystem response. Empirical evidence and the model demonstrate that with multi-factored change, 
the ecosystem response depends on community feedbacks, and that the magnitude of ecosystem response will 
depend on the relationship between plant response to one factor and plant response to another factor. Tradeoffs 
in the ability of plants to respond positively to, or to tolerate, different global change drivers may underlie generalizable 
patterns of covariance in responses to different drivers of change across plant taxa. Mechanistic understanding of 
these patterns will help predict the community feedbacks that determine long-term ecosystem responses. 
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Introduction 

Unpredictability in the functioning of terrestrial ecosys- 
tems underlies large uncertainties in predictions of 
future atmospheric C0 2 content (Friedlingstein et al. 
2006; Denman et al. 2007; Matthews 2007; Dolman 
et al. 2010; Piao et al. 2013). The uncertainty among 
projections of the terrestrial carbon sink over the next 
century is equivalent to 50 years of greenhouse gas 
emissions at current rates (Sitch et al. 2008). Some 
Earth System models incorporate biogeochemical con- 
straints such as those imparted by the nitrogen (N) 
cycle on plant productivity and respiration and thus 
ecosystem carbon exchange (Thornton et al. 2007; 
Sokolov et al. 2008; Bonan and Levis 2010; Bala et al. 
2013), but these models may still overestimate the ability 
of ecosystems to offset human emissions (Piao et al. 
2013). Long-term carbon balance may be sensitive to 
other, slow-acting ecological feedbacks (Friedlingstein 
and Prentice 2010) that affect the magnitude of the 
C0 2 fertilization effect in terrestrial ecosystems. 

Though Earth System models have incorporated more 
sophisticated dynamic global vegetation models to 
represent vegetation shifts, relatively short-term, physio- 
logical mechanisms still underpin plant responses (Smith 
and Dukes 2013). A major source of uncertainty in these 
models is their ability to predict changes in plant commu- 
nities in ways that exert a strong influence on ecosystem 
response to global change (Scheiter et al. 2013), recently 
referred to as the 'community effect' (Polley et al. 2014). 
Plant community shifts in response to global changes 
could initiate key feedbacks that ultimately overwhelm 
short-term physiological responses (Chapin et al. 1997; 
Shaver et al. 2000; Suding et al. 2008; Smith et al. 2009). 
Moreover, the community shifts observed in global 
change experiments may not correspond to those 
expected based on short-term physiochemical responses 
of individual species or functional groups (e.g. Kardol et al. 
2010; Langley and Megonigal 2010; Niu et al. 2010; Dawes 
et al. 2013; Isbell et al. 2013; Smith et al. 2013; Ward et al. 
2013). 

Previous theoretical work has established how changes 
in plant community structure (defined herein as changes 
in species identity, richness or evenness), either from 
shifts of dominance within an ecosystem, immigration 
or extinction, may drive ecological feedbacks on long- 
term time scales (Smith et al. 2009) that are relevant 
for global change predictions. Here, we apply this idea 
to multi-factored change by (i) surveying recent literature 
that highlights the importance of community dynamics in 
long-term, multi-factor field manipulations, (ii) present- 
ing new data from a brackish marsh as a case study, 
and (iii) using a simple community-net primary 



production (NPP) model to illustrate how community 
shifts can mediate ecosystem responses to multiple glo- 
bal changes. We emphasize the marsh as a case study be- 
cause the plant community is simple (richness = 3) and 
dynamic so that it responds rapidly to new conditions. 
The responses observed in the marsh may therefore illus- 
trate mechanisms that apply in other ecosystems that 
have more complex communities such as grasslands, or 
ecosystems that require much more time to equilibrate, 
such as forests. We emphasize the challenge of predicting 
the responses of future ecosystem productivity to ele- 
vated C0 2 in combination with other drivers; however, 
the concepts maybe generally applicable to other combi- 
nations of drivers, as well as other ecosystem responses. 

Community Effects in Response to Single 
Factors 

Single-factor scenarios provide a simple starting point to 
evaluate the importance of the community effect. Rising 
C0 2 should shift communities towards dominance by 
plant species that respond most positively to elevated 
C0 2 (Purves and Pacala 2008; Polley et al. 2012). Typically, 
such a unidirectional community response will have an 
amplifying effect on the ecosystem response and an 
ameliorating effect on the initial resource perturbation. 
For instance, models that allow for vegetation dynamics 
yield stronger negative feedbacks on atmospheric C0 2 
rise than those that have static communities (Purves and 
Pacala 2008). As species that respond strongly to C0 2 in- 
creasingly dominate, enhanced ecosystem C0 2 uptake 
will tend to counteract the resource perturbation, which 
in this case is rising atmospheric C0 2 concentrations. 

Examples from the first generation of global change 
experiments showed in single-factor studies that plant 
community structure responded predictably to one per- 
turbation and had predictable consequences for ecosys- 
tem processes. For instance, N addition strongly shifted 
plant dominance towards nitrophilic species, those with 
a greater capacity to take advantage of extra N, enhan- 
cing total ecosystem productivity in a Minnesota grass- 
land (Tilman 1987). In a brackish marsh, 4 years of 
exposure to elevated C0 2 yielded increasing plant domin- 
ance by C 3 plants over C 4 grasses in mixed plots, enhan- 
cing the overall ecosystem response (Arp et al. 1993). 
These types of findings supported the idea that plant 
community shifts tend to strengthen ecosystem response 
to resource addition. 

However, as global change studies increased in dur- 
ation, encompassing more background environmental 
variability and others included multiple factors, the 
responses became more idiosyncratic. After 18 years of 
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exposure, the same brackish marsh plots showed a weak 
response of plant community to elevated C0 2 compared 
with the community responses to background variability 
in other abiotic factors (Erickson et al. 2007). Similarly, 
elevated C0 2 initially drove a shift in species dominance 
in a New Zealand grassland, but after 11 years, an add- 
itional factor of grazing negated that shift and the ecosys- 
tem C0 2 response (Newton et al. 2014). 

The short- and long-term effects of global change 
drivers may differ for several reasons. For instance, the 
effects on vegetative production may differ from effects 
on reproduction or recruitment (including floral produc- 
tion, seed production, seed quality, seedling establish- 
ment), leading to longer-term demographic shifts that 
are not predictable from short-term growth responses 
(Williams et al. 2007). Moreover, the future will bring 
multiple changes simultaneously, and some factors, like 
N eutrophication, may have stronger effects on plant 
community structure than others, like rising C0 2 (Reich 
2009). What happens to the community feedback when 
multiple factors important to plant physiology and 
growth change at once? 

Community Shifts in Response to Multiple 
Factors 

Relatively few long-term, in situ global change experi- 
ments have manipulated multiple factors (Norby and 
Luo 2004), but several have recently found that the 
combined effect on ecosystems differs substantially 
from the sum of single-factor effects (Table 1). For 
example, a global change experiment in a Californian 
grassland found a C0 2 stimulation of plant growth 
when C0 2 was applied in isolation but not in combination 
with warming or N addition (Shaw et al. 2002; Dukes et al. 
2005), and similarly, C0 2 effects on N 2 0 emissions were 
dampened with other factors (Brown et al. 2012). In a 
hardwood forest, high C0 2 enhanced growth but not in 
combination with ozone (Karnosky et al. 2003) or climatic 
variability (Kubiske et al. 2006). In a tidal marsh, the C0 2 
stimulation of productivity diminished with added 
N (Langley and Megonigal 2010). More generally, the 
magnitude of global change responses decreases with 
an increasing number of drivers (Leuzinger et al. 2011). 

Plant community feedbacks may help explain this 
pattern in ecosystem response. Elevated C0 2 alone 
tends to favour dominance by C0 2 -responsive species 
yielding a more C0 2 -responsive ecosystem (Arp et al. 
1993; Newton et al. 2014; Polley et al. 2014), but addition 
of N can cause unexpected effects mediated by plant 
community shifts (Fig. 1; Langley and Megonigal 2010). 
While elevated C0 2 and N addition tend to have direct 
and positive effects on the productivity of plants in 



isolation, each resource has different effects on plant 
communities by favouring different subsets of species. 
Moreover, N addition generally imparts a much more pro- 
found influence on the community than C0 2 (Reich 2009; 
White et al. 2012). In the tidal marsh, N addition pro- 
moted dominance by C 4 grasses that respond poorly to 
elevated C0 2 , thereby reducing ecosystem responsive- 
ness (Fig. 1, second factor). In years when N had a strong 
effect on dominance by grasses, it also reduced the 
response of NPP to elevated C0 2 (Fig. 2). 

A new generation of decadal-scale global change stud- 
ies is converging on the finding that the strongest effects 
of global change may be mediated by plant community 
shifts, eclipsing other more direct effects on plant physi- 
ology (Kardol et al. 2010; Langley and Megonigal 2010; 
Niu et al. 2010; Yang et al. 2011; Wu et al. 2012; Isbell 
et al. 2013; Newton et al. 2014; Polley et al. 2014). In a 
Tasmanian grassland, plant functional type had a large 
effect on ecosystem carbon cycling (Pendall et al. 2011). 
A recent analysis of N cycle response to warming and 
clipping found that the strongest effects of N on product- 
ivity were driven by species shifts (Niu et al. 2010). In an 
old-field ecosystem, Kardol et al. (2010) found that differ- 
ences in soil processes varied more strongly among plant 
species than between any global change treatments such 
as watering, added N and elevated C0 2 , implying that in 
the long term, even small changes in species abundance 
would hold more sway than the sum of physiological 
changes imparted on individual species. Warming en- 
hanced NPP in an arid grassland until community shifts 
negated the stimulation (Wu et al. 2012). Even though 
N addition strongly increased plant productivity in a 
Minnesota grassland initially, the effect diminished in 
the long term owing to changes in plant diversity (Isbell 
et al. 2013). 

While community shifts may take decades to occur, 
direct physiological responses can occur within seconds 
but may diminish through time because of acclimation 
(Isbell et al. 2013; Smith and Dukes 2013). For instance, 
the short-term physiological response and ecosystem 
productivity response to C0 2 often diminish in the longer 
term (Seileret al. 2009; Lee etal. 2011). Likewise, warming 
of soil causes a short-term increase in respiration that 
wanes over time (Bradford et al. 2008). Therefore, while ex- 
periments are likely to overestimate the direct influence of 
global change, they also may underestimate community 
effects, particularly in plant communities that change 
slowly (Smith et al. 2009). Forest community dynamics 
are particularly difficult to capture owing to long individual 
lifespans. Tree species shifts in response to C0 2 have not 
been adequately examined (Norby and Zak 2011), though 
chronic N pollution is beginning to alter tree species in 
some forests (e.g. Zaccherio and Finzi 2007). 



AoB PLANTS www.aobplants.oxfordjournals.org 



© The Authors 2014 



3 



Table 1. Long-term, multifactor, C0 2 studies that report on plant community feedbacks. 



Site 



C0 2 treatment Other Duration C0 2 effects on NPP 

treatments (years) 



MN grassland Ambient + 180 ppm 
BIOCON 



Swiss pasture Ambient + 240 ppm N 



MD Marsh 720 ppm 



13 



10 



CA grassland 680 ppm 
Jasper Ridge 



Tasmanian 550 ppm 
grassland 

WY Prairie 600 ppm 



Danish 510 ppm 

heathland 



Temp, N, H 2 0 4 



Temp 



Temp 



Temp, H 2 0 3 



13 % (R&H, 2012) 



Increased 
belowground but 
not AGB. Decreased 
N concentration 

10 % over 4 years 



No significantly effect 
onshoot or root 
production 



Plant community response 
to C0 2 only and/or effects on 
response 

C0 2 increased biodiversity, 
yielding greater C0 2 
stimulation 



More T. repens at elevated 
C0 2 (33 % instead of 21%) 
at expense of L perrene 
(Hartwig et al. 2000) 

No 



Forb abundance was lower 
(Zavaleta et al. 2003) 



No 



Multiple-factor community 
feedbacks 

N addition enhanced C0 2 
response of ANPP. Interactive 
effects independent of 
community richness or 
composition 

High N addition enhanced C0 2 
response of AGB 



60 % increase in ANPP C0 2 increased C 4 dominance 
over 4 years from 25 to 39 % 

Peak total AGB increased Favoured C 3 grasses (34 % 
an average of 33 % greater AGB than control) 
(but not during wet over C 4 (28 % greater) 
year). 

Increased biomass of 
D. flexuosa (grass) in 1 of 3 
years. No significant effect 
for C. vulgaris (evergreen 
dwarf shrub) 



N addition strongly shifted plant 
community towards G, 
dominance, reversing C0 2 
effect 

Adding warming or precipitation 
to C0 2 treatment increased forb 
abundance, while adding N 
decreased it, no C0 2 effect 
with N 

Warming decreased C 4 
dominance. No interaction 
reported 

Warming favoured C 4 grasses 
reducing C0 2 effect 



Citation 



Isbell etal. (2013), 
Reich and Hobbie 
(2012) 



Schneider et al. 
(2004), Hartwig 
et al. (2000) 

Langley and 
Megonigal (2010) 



Less D. flexuosa biomass in 
warmed + drought plot 
compared with drought plots 
but C0 2 counterbalanced this 
decrease in three-factor plots 



Dukes et al. (2005), 
Zavaleta et al. 
(2003) 



Pendall et al. (2011) 



Morgan et al. (2011) 



Kongstad et al. 
(2012) 
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Figure 1. The pathways of influence by which elevated C0 2 and N 
may drive an ecosystem function. Symbols represent the effect of 
added resources on each relationship. The black arrows represent 
short-term, physiological effects of N addition on the C0 2 fertiliza- 
tion effect, which are generally found to be positive for individual 
plants. The red arrows represent the effect of N-driven community 
shifts, which may counteract or enhance the physiological effects 
depending on how they covary with the other pathways of influence. 
Elevated C0 2 may also affect community change, but the effects are 
typically not as pronounced as with added N. 
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Figure 2. Annual NPP response to elevated C0 2 (elevated NPP- 
ambient NPP) by annual grass contribution to community domin- 
ance (percent grass biomass in elevated C0 2 treatment or elevated 
C0 2 + N treatment for N-fertilized plots) from a tidal wetland. Data 
are grouped by N treatment (no N: white circles, added N: black 
circles) from 2006 to 2012. Methodological details are as described 
in Langley and Megonigal (2010); data presented here show add- 
itional years of the study (2010-12). 

Thus, the key to predicting long-term ecosystem 
responses and responses to multiple factors may lie 
in predicting long-term plant community changes. For 
instance, while a diminished ecosystem response to 
multiple factors may be partially explained by different 
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treatments counteracting each other (i.e. warming may 
dry soil while water addition will restore moisture) (Luo 
et al. 2008; Leuzinger et a/. 2011), mounting evidence 
indicates that in the long term, changes in plant commu- 
nity structure could ultimately govern ecosystem re- 
sponses (Figs 1 and 2; Smith et al. 2009; Langley and 
Megonigal 2010; Niu et al. 2010). Yet, accurately predict- 
ing these phenomena requires understanding commu- 
nity responses mechanistically. 

Tradeoffs in Plant Responsiveness 
to Different Perturbations 

A cornerstone of modern ecological theory is that plants 
exhibit tradeoffs in resource strategy (Tilman 1988; Craine 
2009), such that optimization for acquisition and use of 
one resource precludes optimal acquisition and use of 
another. A simple example of a tradeoff is that of plant 
allocation among plant organs. Energetic limitations 
dictate that the same resources cannot be allocated 
simultaneously to multiple organs, while phylogenetic 
limitations appear to dictate that a single plant has lim- 
ited plasticity in shifting allocation— that is, plant species 
are hardwired to a narrow range of allocation strategies 
(Reynolds and D'Antonio 1996; Muller et al. 2000; Craine 
and Dybzinski 2013). For example, N addition induced 
only small changes in plant allocation patterns compared 
with the inherent variability across 27 herbaceous species 
(Muller et al. 2000). 

These differences in allocation strategy could help 
explain how plants respond to global changes, particularly 
those that alter resource availability. For instance, 
plants that exhibit low root-to-shoot biomass ratios are 
optimized for relatively high nutrient conditions, and 
should be better positioned to respond positively to soil 
resource addition than plants that have high root- 
to-shoot ratios. Other tradeoffs in physiology or allocation 
may play a role in resource strategy. For instance, a trade- 
off may exist in plant metabolic pathways such that a 
positive response to elevated C0 2 may interrupt the 
efficient incorporation of nitrate-N into a usable organic 
form (Bloom et al. 2010). If these tradeoffs result in 
a predictable relationship between the magnitude of 
plant growth response to elevated C0 2 and that to 
added N, then this knowledge could help constrain our 
predictions of ecosystem responses to global change. 

Evidence for a tradeoff between C0 2 and N responses 
exists across functional groups. For instance, legumes 
respond most strongly to C0 2 compared with other func- 
tional groups (Jablonski et al. 2002) and are among 
the weakest responders to N addition (Xia and Wan 
2008). At the other end of the spectrum, grasses tend 
to exhibit a very weak C0 2 response (Jablonski et al. 



2002; Ainsworth and Long 2005), but respond strongly 
to N enhancement in many ecosystems (Pennings et al. 
2002; Langley and Megonigal 2010). Slow-growing plants 
generally respond more positively than fast-growing ones 
to elevated C0 2 (Ali et al. 2013), while the opposite is true 
for responses to N addition. Still, enormous variability ex- 
ists within functional groups, making it difficult to 
generate useful functional group parameterizations for 
models. Again, these tradeoffs in response may occur 
among other resources or conditions as well. Recent 
evidence, for example, suggests differential responses 
to C0 2 and water, and the same responses to C0 2 and 
warming (Dijkstra et al. 2010; Dieleman et al. 2012). 

A Simple Plant Community Model 

To examine the potential importance of plant community 
shifts for ecosystem response we simulated a simplified 
ecosystem with 30 species of initially equal biomass. We 
assigned realistic, randomized distributions of biomass 
responses to each plant species modelled after observed 
distributions of responses. Carbon dioxide response ratios 
(R = treated/control) across species are generally nor- 
mally distributed with a mean of 1.25 and a standard 
deviation of 0.35 (Poorter and Navas 2003). Nitrogen sti- 
mulations are typically log-normally distributed (Xia and 
Wan 2008) with a mean of 1.29 and a standard deviation 
of 0.37 (LeBauer and Treseder 2008). (Note that model 
outputs are qualitatively insensitive to the effect sizes 
used.) We assigned a random value (P) from 0 to 1 to 
each of 30 species. We used Pas a cumulative probability 
to calculate the corresponding C0 2 response according to 
a normal distribution (/x = 1.25, a= 0.35). 

Then, we generated combined (C0 2 + N) responses 
for three types of relationships between C0 2 and N 
responses: negative, positive and independent, repre- 
senting three ways that responses to each driver may 
covary. For 'negative', the randomly generated probability 
was inverted (1-P) to generate the N response of 
each species using a lognormal distribution (jx = 1.29, 
a= 0.37). Therefore, a species with a large C0 2 response 
would have a relatively small N response and vice versa. 
The negative relationship represents a tradeoff wherein 
plants that respond strongly to C0 2 do not respond 
strongly to N. For 'positive', each species probability (P) 
was converted to an effect size for N using the lognormal 
distribution of N responses (/x= 1.29, a= 0.37) such that 
a species that responds most strongly to C0 2 also re- 
sponds most strongly to N. A positive relationship could 
occur if the capacity to respond to both N and C0 2 was 
related to some other factor. Some have suggested, for 
instance, that plants with high relative growth rates 
should be more suited to take advantage of increased 
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resource availability, though this has not been observed 
for elevated C0 2 (Poorter and Navas 2003). The 'inde- 
pendent' scenario assumes that the N response of a 
species is unrelated to the C0 2 response, so random 
probabilities were generated independently for the N 
responses and then converted using the N response dis- 
tribution (/jl= 1.29, cr= 0.37). The response ratios were 
maintained for each species throughout the simulation. 

We calculated a species biomass response ratio (R) for 
each species in each scenario. For single-factor scenarios 
(C0 2 only and N only), we used R for each species based 
only on either C0 2 or N responses. For the C0 2 + N scen- 
arios (independent, positive and negative) we calculated 
additive and multiplicative effects of C0 2 and N. Additive 
was estimated as R x = Rqo 2 + ^Nitrogen- Multiplicative 
stimulation was calculated as R x = Rco 2 x ^Nitrogen- 
These two methods represent two ways C0 2 and N may 
interact in ecosystems (Reich et a/. 2006). Because the 
additive and multiplicative simulations yielded qualita- 
tively similar responses, we limit further discussion to 
the multiplicative. We weighted each plant response 
according to its relative abundance (A x ) to calculate 
total ecosystem biomass response (ER): 

ER = £(A X R X ) 

Therefore, the first-generation ecosystem response is 
equivalent to the average of individual plant responses 
because the distribution of species was initially even. 
For subsequent generations, we allowed species relative 
abundance in a given generation (t) to shift according to 
its relative biomass response in the previous generation 
(t - 1): 

Axm = Rx(t-i)/ER(t-i) 

Here we assumed that community shifts follow biomass 
responses. This assumption may not hold as, for 
instance, reproductive response, which may influence 
long-term demography, can differ from the growth 
response (Williams et a/. 2007; Way et al. 2010). Still, 
there is evidence that long-term plant community shifts 
may respond consistently over time to some drivers like N 
eutrophication more strongly than to others like elevated 
C0 2 (Reich 2009; Langley and Megonigal 2010). To ac- 
count for the possibility that community shifts follow N 
responses but not the combined C0 2 x N responses, we 
conducted a simulation using just the single-factor R x 
representing the N response. 

Because no species are added and none go extinct, 
the model represents physiological responses and shifts 
within the community, or species reordering (Smith 
et al. 2009), but not exchange with external species 



pools. However, species relative abundance is allowed 
to drop infinitesimally close to zero, so that species may 
effectively go extinct in some cases. 

To constrain the upper limit of ecosystem response we 
imposed a light limitation according to an established 
relationship between leaf area index (LAI) and canopy 
assimilation given a high light environment given the 
following equation (Bonan 2008): 

"canopy — l/^max /K) x In [(A m ax/e + F 0 )/(A max /e 

where A canopy is the ecosystem level C0 2 assimilation per 
ground area used as a proxy for plant growth; A max is the 
light saturated leaf-level assimilation (20 |j,mol nrT 2 s~ a ); 
e is the light use efficiency (0.06); K is the light extinction 
coefficient (0.5); F 0 is the solar radiation at the top of the 
canopy (1000 ixmol itT 2 s~ a ); and L is the LAI, which was 
set initially at 2.5 and allowed to vary through time. Then, 
we allowed shifts to occur for 14 additional generations to 
look at the trajectory of ecosystem response with shifting 
community structure (Fig. 3). We repeated the simulation 
with re-randomized response distributions 20 times to esti- 
mate mean response and error. 

Though this exercise makes many simplifying assump- 
tions, the results illustrate some important points. If the 
ultimate community shifts follow initial, short-term 
responses, then the community feedback effect will 
amplify individual physiological responses to an en- 
hanced resource in isolation to other global change fac- 
tors or other disturbances ('C0 2 only' and 'N only' 
scenarios; Fig. 3, top panel). The ultimate ecosystem 
NPP response roughly doubles the initial stimulation, 
intuitively approaching the stimulation of the strongest 
responding species. This finding agrees with initial results 
from many ecosystem studies described above. 

With multiple interacting factors, the magnitude of 
the initial ecosystem response (at Generation 1) and es- 
pecially the community effect (apparent in subsequent 
generations) depends on how plant responses to one fac- 
tor relate to that of another factor. The existence of 
a negative relationship between C0 2 response and N 
response across plant species, which would arise if trade- 
offs exist in the ability of plants to respond to enhance- 
ment of each resource, will reduce total ecosystem 
productivity in the long term below what would be pre- 
dicted if C0 2 and N responses are independent or posi- 
tively related. Interestingly, the model indicates that 
with negative covariance, community feedbacks in 
response to C0 2 and N in combination could lead to an 
ecosystem productivity stimulation that is smaller than 
with C0 2 only (Fig. 3). The short-term responses do not 
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Generations 

Figure 3. Ecosystem productivity response (top panel) and relative 
difference in productivity (bottom panel) for simulated plant com- 
munities for single factors (C0 2 and N) and combined responses 
(C0 2 + N) according to the relationships between C0 2 response 
and N response across species (independent, positive or negative). 
Ecosystem productivity is expressed as a factor of the control with 
ambient C0 2 and no added N, which is set to 1. The first-generation 
response (highlighted in grey box) represents the physiological effect 
while subsequent generations represent the influence of dynamic 
communities. For instance, while C0 2 alone here yields a response 
factor of 1.24 in the first year, ensuing community shifts favouring 
C0 2 -responsive species could increase the response to 1.49 by year 
2015. 



depend as strongly on the type of covariance (from 42 to 
56 % stimulation, Generation 1; Fig. 3) as the long-term 
responses (from 46 to 86 % stimulation, Generation 15; 
Fig. 3). 

We have made the point above, however, that long- 
term community shifts do not always follow initial physio- 
logical responses. If the ultimate community shifts occur 
randomly and independently of the physiological stimu- 
lation, then the ultimate stimulations would resemble 
the initial stimulation (Generation 1), and community 
shifts would be unimportant for predicting long-term 
ecosystem response. It is more likely that community 
shifts are related to physiological responses to global 
change drivers, though the shifts may not occur in a 
manner as simple as represented in this model. For in- 
stance, N addition tends to have a stronger influence on 
community structure than does elevated C0 2 (Reich 
2009; Langley and Megonigal 2010). Therefore, we may 
expect that plant communities exposed to both C0 2 
and N would follow the initial N response more than the 



initial C0 2 response. When community shifts followed N 
response, the results were dampened but qualitatively 
similar to when the community followed the combined 
C0 2 x N response (results not shown). 

The model indicates that failing to properly account for 
plant community shifts could engender large errors in 
estimates of long-term ecosystem productivity. In the 
case of multi-factored global change, the community 
feedback effect depends strongly on the type of covari- 
ance plant responses exhibit among different factors. 
For instance, failing to account for negative covariance 
between C0 2 and N response could lead to a drastic 
overestimate in ecosystem response. Recognizing these 
patterns could help improve predictability of long-term 
plant responses. Indeed several examples from the 
recent literature point to such shifts in community struc- 
ture driven by one factor diminishing the ecosystem 
response to another (Table 1) indicating that capturing 
these types of interactions in models will help predict 
future response of ecosystems to global change. 

Predicting the Global C0 2 Fertilization 
Effect 

A decade ago, global climate models systematically 
overestimated the ameliorating effects of terrestrial 
ecosystems on the future carbon cycle by excluding the 
nutrient limitation of the C0 2 fertilization effect. Several 
more recent models have now incorporated this stoichio- 
metric limitation (Sokolov et a/. 2008; Thornton et a/. 
2009; Bala et al. 2013). Interestingly though, as models 
refine representation of the N cycle, some show that 
N may become less limiting in the future, owing to accel- 
erated mineralization in response to warming and projec- 
tions of increasing N deposition. Accounting for enhanced 
N availability, therefore, partially restores the large C0 2 
fertilization effect (Sokolov et al. 2008; Thornton et al. 
2009; Zaehle and Dalmonech 2011). Evidence from long- 
term C0 2 studies suggests that N scarcity limits the C0 2 
stimulation of productivity in temperate ecosystems 
(Norby et al. 2010; Reich and Hobbie 2012). However, N 
addition exerts a strong influence on plant diversity 
and species composition (Tilman 1987; Isbell et al. 
2013) that must be considered when extrapolating to a 
future with widespread and heterogeneous increase in 
N availability. 

It has been suggested that Earth System models, many 
of which represent biodiversity coarsely, overestimate 
the negative effects of global change on ecosystem pro- 
cesses because they exclude the mitigating effects of 
plant species shifts (Purves and Pacala 2008) apparent 
in diversity experiments (Tilman et al. 2006). For instance, 
higher plant biodiversity diminished the negative effects 
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of drought on plant productivity (Tilman etal. 2006; Reich 
2009). However, in predicting an enhanced terrestrial 
carbon sink in the future (Pinsonneault et al. 2011; Bala 
et al. 2013), the expectation is reversed — rather than rely- 
ing on stable ecosystem functioning, society stands to 
benefit from a dynamic ecosystem function, specifically, 
increasing net carbon uptake through time. If this dra- 
matic change in ecosystem behaviour does not persist, 
then the atmospheric C0 2 concentration will rise more 
than expected, and climate change will be more severe 
than presently predicted. 

Conclusions 

We have focused on how the interactive effects of two 
added resources may affect plant communities to drive 
ecosystem function. Yet, the recognition of community 
shifts in global change response applies to other factors 
besides enhanced resource availabilities (De Marco et al. 
2013). For instance, if ozone tolerance were negatively 
related to C0 2 responsiveness across plant taxa, then 
ozone exposure will favour species that exhibit a slighter 
C0 2 response, thereby dampening the ecosystem-level 
C0 2 response. The combination of C0 2 and N could like- 
wise favour a unique subset of species that can take 
advantage of the combination of high C0 2 and high 
N. For example, some invasive species have been shown 
to respond more strongly to both elevated C0 2 and 
N enrichment than non-invasives (Dukes et al. 2011; 
Mozdzer and Megonigal 2012). This general propensity 
for invasives to take advantage of high resource availabil- 
ity (Davidson etal. 2011) may indicate that they are not as 
constrained by the same tradeoffs as other plants. The 
degree to which community will shift and ecosystem 
productivity will respond to future change will partly 
depend on the degree to which individual plant responses 
to global change drivers covary across taxa. 
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